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Abstract 
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1. Introduction 

Measuring purely leptonic decays of heavy mesons provides the most clear way for the 
determination of weak decay constants of heavy mesons, which connect the measured quantities, 
such as the BB mixing ratio, to CKM matrix elements Vcb, V^h- However, currently it is not 
possible to determine /b, fs^ fo^ and /d* experimentally from leptonic B and Dg decays. For 
instance, the decay rate for is given by |l| 

Because of helicity suppression, the electron mode —>■ e^u has a very small rate. The relative 
widths are 10 : 1 : 2 x 10~^ for t^v^^'^v and e^v final states, respectively. Unfortunately the 
mode with the largest branching fraction, r+i/, has at least two neutrinos in the final state and 
is difficult to detect in experiment. So theoretical calculations for decay constant have to be 
used. The factorization ansatz for nonleptonic decay modes provides us a good approximate 
method to obtain nonperturbative quantities such as form factors and decay constants which 
are hardly accessible in any other way 0, |^ . 

There are many ways that the quarks produced in a nonleptonic weak decay can arrange 
themselves into hadrons. The final state is linked to the initial state by complicated trees 
of gluon and quark interactions, pair production, and loops. These make the theoretical de- 
scription of nonleptonic decays difficult. However, since the products of a 5 meson decay are 
quite energetic, it is possible that the complicated QCD interactions are less important and 
that the two quark pairs of the currents in the weak Hamiltonian group individually into the 
final state mesons without further exchanges of gluons. The color transparency argument sug- 
gests that a quark-antiquark pair remains at a state of small size with a correspondingly small 
chromomagnetic moment until it is far from the other decay products. 

Color transparency is the basis for the factorization hypothesis, in which amplitudes fac- 
torize into products of two current matrix elements. This ansatz is widely used in heavy quark 
physics, as it is almost the only way to treat hadronic decays. 

In this paper we consider the way how to determine weak decay constants f^^ and f^* 
under factorization ansatz including penguin effects. In section 2 we discuss the way how to 
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extract the unknown parameter \VcbFf^{0) \ from the branching ratio of the semileptonic decay 
—>■ D~^iu. In order to check the vahdity of the factorization assumption, we study the 
nonleptonic two body decays, B Dp, Dn and DK^*^ in section 3. In section 4 we calculate 
fns and fn* from B^ D^Dj^*^ decay modes. In our analysis we improve the previous 
analysis by considering the QCD-penguin and Electroweak-penguin effects of about 13 % 
for B DDs and 4 % for B DD*, which are not negligible as discusssed in Also we 
follow the gauge independent approach to calculate the effective Wilson coefficients which was 
studied by pertubative QCD factorization theorem 0. 



2. Semileptonic Decay D^l u 

From Lorentz invariance one finds the decomposition of the hadronic matrix element in 
terms of hadronic form factors: 



<D+{Pd)\J,\B\pb) > 



{Pb+Pd] 



2 2 



m% — m 



D 



(2) 



where = 07^6 and = {pb — Pd)/x- In the rest frame of the decay products, Fi{q^) and 
Fo(g^) correspond to 1~ and 0+ exchanges, respectively. At = we have the constraint 



Ff ^(0) = F(f"(0) 



(3) 



since the hadronic matrix element in (0) is nonsingular at this kinematic point. 

The distribution in the semileptonic decay B^ D^l^v is written in terms of the 
hadronic form factor F^^{q^) as 

c/r(50 ^ D+l-u) 



V^,\'[K{q')f\Fr{q'''' 



where the q^ dependent momentum K{q'^) is given by 
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K{q' 



2m I 



ttib + — Q ) — AmBiTLj^ 
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(4) 
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In the zero lepton mass limit, 0<g^< {ttib — moY 



For the dependence of the form factors, Wirbel et al. [0] assumed a simple pole formula 
for both -Fi(g^) and Fo^q^) (pole/pole): 



Fi(g^)=Fi(0)/(l 



q^ 



m 



2 >■> 

Fi 



W)=i^o(0)/(l-^) 



m 



2 

Fo 
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with the pole masses 



mp, = 6.34 GeV, = 6.80 GeV. (7) 

Korner and Schuler ^ also adopted the same dependence of Fi{q^) and F^i^q^) given by (||) 
and d^). On the other hand, the heavy quark effective theory (HQET) gives in the m^ c oo 
limit the relation between Fi{q^) and -fo(?^) given by |l^ 

-,2 



Fo{q' 



1 - 



Fi{q' 
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The combination of (^) and suggests that F^lq"^) is approximately constant when we keep 
the simple pole dependence for Fi(g^). Therefore, in this paper, as well as the above (pole/pole) 
form factors, we will also consider the following ones (pole/const.): 
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Fi(g^)=Fi(0)/(l 
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2 ^' 
Fi 



Fo(g^)=Fo(0), 



(9) 



with 



mp, = 6.34 GeV. (10) 

By introducing the variable x = q'^/m'^, which has the range of < x < (1 — in the 

zero lepton mass limit, is written as 



D+l-y) 



dx 
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Then the branching ratio B{B^ — > D^l v) is given by 



^ D+l-v) = ( 



GplTil 2 ITT-B 2 



v/2 ' 1927r2 
2.221 X 10^ \VcbFf^ 



v^bFrm'^i 



X /, 



(12) 
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where the dimensionless integral I is given by 



(1- 



dx 



1 H — r — a; r — 4— f- 



1 3/2 



(1 



0.121 



(13) 



In obtaining the numerical values in (|T2D and ([T3|), we used the following experimental results 
|ni : mz) = = 1.869 GeV, tub = trbo = 5.279 GeV, Tb = Tbo = 4.219 x 10"^=^ GeV 

(tbo = (1.56 ±0.06) X 10-12 g^j^j ^ 1 igg 39(^2) x 10"^ GeV~2. Since B{B^ D+l-y) = 
(1.78 ± 0.20 ± 0.24) X 10"^ was obtained experimentally, the value of |Kfe^f^(0)| can be 
extracted from (|12|). Following this procedure, we obtain 



VcbFl^^{0) \ = (2.57 ±0.14 ±0.17) x 10" 



(14) 



In the calculations of the next sections, we will use \ VcbFf^{0)\ = (2.57 ± 0.22) x 10 ^ which 
is given by combining the statistical and systematic errors in (|T^) . 



3. Test of Factorization with D^p and 5 D^tt , and Prediction 

of Branching Ratio B{B^ D+R-^*^) 

In general the test of factorization, independent of the numerical values of ai, a2 and of the 
CKM parameters |V"cfe| or \Vub\, can be carried out by considering the ratios of rates for two 
Class I or Class II i?-meson hadronic two-body decays. On the other hand, we can also use the 
relation between the semi-leptonic decays and the non-leptonic decays with ai and 02 given by 
other sources. In our analysis we use the latter one. 

Let us start by recalling the relevant effective weak Hamiltonian: 

Kff = ^V,bV:,[C,{^^)0, + CMO2] + H.C., (15) 

Oi = {dVPu){cVpb), 02 = {cVfu){dVpb), (16) 

where is the Fermi coupling constant, Vd, and Vud are corresponding Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements and Fp = 7p(l— 75). The Wilson coefficients Ci(/i) and C2(/i) 
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incorporate the short-distance effects arising from the renormahzation of TieS from fi = mw 
to fi = 0{miy). By using the Fierz transformation under which V — A currents remain V — A 
currents, we get the following equivalent forms: 



(C2 + Y^i)02 + Ci(cr^T"n)(dr,T'^6), (17) 



where iVc = 3 is the number of colors and T"'s are S'f/(3) color generators. The second terms 
in ([TtD involve color-octet currents. In the factorization assumption, these terms are neglected 
and Tioff is rewritten in terms of "factorized hadron operators" [0: 

H,^ = ^V,,V:^[a^[dTPu]H[(T,h]H + a2[cTP^^^^ + H.C., (18) 

where the subscript H stands for hadronic implying that the Dirac bilinears inside the brackets 
be treated as interpolating fields for the mesons and no further Fierz-reordering need be done. 
The phenomenological parameters ai and 02 are related to Ci and C2 by 

ai = Ci + ^C2, a2 = C2 + ^Ci. (19) 

From the analyses of A.J. Buras [0, the parameters ai and 02 are determined at NLO calcu- 
lation in the NDR scheme as 

ai = 1.02 ±0.01, a2 = 0.20 ±0.05. (20) 
For the two body decay, in the rest frame of initial meson the differential decay rate is given 

by 

I [(M^ - (mi ± m2)^)(M^ - jm, - m^)')]'/' 

iPil 2M ' ^ ^ 

where M is the mass of initial meson, and mi (m2) and pi are the mass and momentum of one 
of final mesons. By using (|^), (|18D and < 0|r^|p(g, £:) >= e^{q)mpfp, (|2T|) gives the following 
formula for the branching ratio of B'^ D^p^: 

B{B p) = (^-^j IK.I i^T7«i;;^l^^^^i K)l 



X 



1 - 



rriB 



1 - 



rriB 



3/2 



13.25 X K)r X 



1.02 



(23) 



In obtaining the numerical values in (p3D, we used the experimental results given below (JT 



766.9 MeV, L = /„+ = 216 MeV, and Kd = 0.9751 For the value of Oi we 



used the value given in (|20|) . Then, by using the formula ( p3|) with the values of \VcbFQ^{0)\'^ 
(F(f^(0) = Ff^{0)) given in ([l^) , we obtain the branching ratio B{B^ D^p^) presented in 
Table |l|. 

For the process — > D^K*^, by using < 0\T^\K*{q,e) >= e^{q)mK* fx*, we have 



B{B^ D+K* 



( GpniW 2 1 TUB 2 fh 
V V2 / 167r Vb m% 



m 



K* 



1 - 



,m£) — rriK* ^2 



TUB 
?BD/^2 m2 



0.67X |V;,Ff^(m^Or X 



ai 



niB 
2 



3/2 



1.02 



(24) 



where we used tuk* = niK*- = 891.59 MeV, fx* = fx*- = 218 MeV, and Ks = 0.2215 [|Tl[. By 
using (p^ ) with \ Vcb Ff^(0)| ^ in (p!^), we obtain the branching ratio B{B^ — > D^K* ) presented 
in Table |I[ 

By using (^), (|T8|) and < 0|r^|7r(g) >= iq^f-K, ( PI) gives the following formula for the 
branching ratio of the process B^ —>■ D^n^: 



B{B^ 



D + 71- 



(GfTUb] I |2 1 mB 2 L IT. pSD/ 2n|2 

V V2 / 167r Fb m| 



x(l 



m^x2 



1 



)^)(i-( 



niBi — 
ttlb 



f) 



1/2 



5.42X |K6Fo''''K)l'x 



ai 
1)2 



(25) 



where we used = mT^~ = 139.57 MeV and = f^^- = 131.74 MeV |[TT|. By using the formula 
(H) with the values of |V;bF(f^(0)p (F(f^(0) = Ff^{0)) in (|T|), we obtain the branching ratio 
B^ D^n^ presented in Table |1|. 

For the process B^ D^K^, by using < 0\Tf^\K~{q) >= iq^fx-, we have 

c2 



B{B^ 



(^) 



2 1 "^-B ^2 iTA t:^BD/2 \\2 



167r Fs ^ TOb 



[m 



Kj 



X 1- 



m 



D 



m 



B 



1 - 



.rriD + mK.2 



rriB 



TUB 



1/2 



0.41 X |K,Fo^^(m^)|^ X 



(26) 



where we used = tuk- = 493.68 MeV, fx = fK+ = 160.6MeV ITTl. By using (|2|) with 
^{0)1"^ in (14), we obtain the branching ratio B{B^ — > D^K^) presented in Table |l[ 



It seems that the factorization method works well in 
branching ratios : 



D^Ti , p decays. We predict 



S(S° 
i3(5° 



2.7- 10"^ ■ 
i 4.6 ■ 10"^ 



ai 

L02 

ai 

L02 



(27) 



which is certainly reachable in near future. 



4. Determination of Jd* and fo from D^D^ * and D^D^ 



The effective Hamiltonian for Ai? = 1 transitions is given by 

i/eff = %[Vu,V:^{C,0'^ + C^O") + V,,V:^{C,Ol + C,0^) - Vt,V:^ C^Ol (28) 

where q = d,s and Cj are the Wilson coefficients evaluated at the renormalization scale /i, and 
the current-current operators Oi''^ and 02^'^ are 

O" = {Uaba)v-A{qf3U(3)v-A = iCaba)v~A{q(3Cf3)v^A 

O2 = {upha)v-A{qaUp)v-A = {cpba)v-A{qaCp)v-A, (29) 

and the QCD penguin operators O3 — Og are 

O3 = {qMv-A J2Wf^)v-A O4 = {qph^)v-A J2i^Mv~A 

q' q' 

O5 = {qaK)v-Aj2(^Wf^)v+A Oq = {qpb^)v-Aj2(^Mv+A- (30) 

q' q' 
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The electroweak penguin operators O7 — Oiq are given by : 

3 3 

q' q' 

3 3 

Oq = iqaba)v-AY.2^<l'Wfs)v-A Og = (gA)y-A ^ 2 (^a^^) V-^- (31) 

<?' 9' 
In (^) we consider the effects of the electroweak penguin operators, however, we neglect the 
contribution of the dipole operators, since its contribution is not important in this work. 

When we take rrit = 174 GeV, rrib = 5.0 GeV, a,{M,) = 0.118 and aem(M^) = 1/128, the 
numerical values of the renormalization scheme independent Wilson coefficients at /i = nib 
are given by [|I^ 

Ci = -0.3125, C2 = 1.1502, 

C3 = 0.0174, = -0.0373, C5 = 0.0104, Cq = -0.0459, 
C7 = -1.050 X 10-^ Cs = 3.839 x 10-^ 

Cg = -0.0101, Cio = 1.959 X 10-\ (32) 
The effective Hamiltonian in (28) for the decays D^D~^*^ can be rewritten as 
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Heff = [K^K: (Ct^O'i + CI" 01) - V,bV,: Y.Cf^O,i (33) 



where are given by ||T5 



with 



C-i — C-i, (_/2 — O2, O3 — C-3 — i^/iVc, U4 — C-4 + -Ts, 

_ /=y p / AT r"^ff — H j_ p r"^ff — n \ p r<f^ff _ 

C-5 — C-5 — i^/iVc, (_/g — C-6 + -Ts, 07 — U7 + -Tg, Ug — Og, 

Cg'^'^ = C*9 + Pe, CiQ^ = ClQ. (34) 



Ps = - G{mg,q^,fi)]C2in), 

on 9 

A = ^[^-G(m„g^/x)](3C'l(/.) + C2(/i)), (35) 

G(m„, q , jj.) = —4 / a;( 1 — x) m( — jax, 

Jo 



(36) 



where q denotes the momentum of the virtual gluons appearing in the QCD time-hke matrix 
elements, and A^c is the number of colors. Assuming = itlI/2, we obtain the analytic formular 
for GijUq, g^, /i) : 



rUb 



3 V8 



+ 



10 

9" 



(2 + y)v/r^ 



In 




+ in 









(37) 



with y = Sniq/ml. 

By considering the non-factorizable term contributions, the relation between the effective 
coefficients a^-^-^ and the Wilson coefficients in the effective Hamiltonian are given by 



eff 
hi 



r^ff 



+ 



1 



.a 



eff 



ff^2i-l, 



eff _ (^eff I <^ 



eff 



at; 



(38) 



where i = 1, 



5, and the non-factorizable effects are absorbed into the N^-^'^ by 
1 1 



- AT + ^" 



(39) 



In order to simplify the notation, we will use the notation a* instead of a^^^ in the below. 

In usual factorization approach, when we consider the off-shell momentum of the external 
quark line, the effective Wilson coefficients has the ambiguities of the infrared cutoff and gauge 



dependence. As stressed by |[16[, the gauge and infrared dependence always appears as long as 
the matrix elements of operators are calculated between quark states. Recently this problem 
was sloved by pertubative QCD factorization theorm by using the on-shell external quark. 
By following their approach and inserting the values for rriq = mc{^) = 0.95 GeV, we get the 



values C' 



eff/ 



1 ~ 10) for 6 ^ c given in Table For different combinations of N^^^ = 2, 3, 



and 5, the values of the effective coefficients aj(2 = 1 ~ 10) are shown in Table ^. Here 
{Nc)ll,lr = 3 corresponds to the naive factorization approximation without considering non- 
factorizable contributions. 

The decay amplitude A{B'^ D^D^) =< D^D^\Hc{{\B'^ > is given as follows: 

Gp. . . m'' 



V2 



[VcbV*^ai - VtbVtKai + aio + 2(a6 + ag) 



-)] 



~ -j=VcbV*^ ai RdDs Ma 



(40) 
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where 



R 



DDs 



m 



Ds 



^ _^ (Q4 + aio) _^ ^ (qg + as) 

ai ai (mb - mc)(m^ + 



and 



Ma =< D;\sYl5c\0 >< D+\c-f^b\B' >= -ifuAml - ml)F^''{ml^ 
On the other hand, we have 

G p . 



(41) 



(42) 



V2 



[V;6K>1- W.(«4 + «10)] Mb 



where 



and 



R 



G 

-^VchV*^ ai Rdds* Mb 



(04 + aio) 



DDs* 



ai 



(43) 



(44) 



Mb =< >< D^\cj,b\B' >= mn*fDMl) " (pb + Pd)]^^ ""(m?,.). (45) 

We can estimate the penguin contributions for each process, for exapmle, in the case of 
Nll = 2 and Nlr = 5 : 



For ^ D+D: 



For D+D;-; 



At 



A 



■T 



(04 + Qio) _^ 2 '^"'^ "'^^ 



m 



Ds 



A, 



An 



ai 



(04 + flio) 



ai {mb-mc){mc + ms) 



ai 



3.9% 



13.1%(46) 



(47) 



where At{Ap) stands for the amphtude of tree diagram (penguin diagram). Here we used the 
values mc{mb) = 0.95 GeV and ms{mb) = 90 MeV. Therefore, the penguin contributions affect 
the extraction of the decay constants fo^ and fp* . The penguin contributions for B DDs is 
more than three times of those for B DD^. 

From (^0|) and (|43|) the decay constants are given by 



= (0.87 X 10"^ GeV) • 



/^^ = (0.64 X 10"^ GeV) 



B{B^^D+D-*) fim 



\V,bF,^^{ml 



ai 



R 



DDs* 



B{B^^D+D-) /1.02 



\VcFr{ml^)\ 
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R 



(48) 



DDs 



Browder et al. ||T^ presented the following experimental results for the branching ratios: 



B{B^ D+D-) 



1.14 ±0.42 ±0.28) X 10^ 



;i.l4±0.50) X 10 



-2 



(0.74 ± 0.22 ± 0.18) X 10"^ = (0.74 ± 0.28) x 10 



1-2 



(49) 



where we combined the statistical and systematic errors. From (|4l|) , (^) , (^) , and (^91), we 
obtain the results which are obtained by including the penguin contributions: 

/d* = 346 ± 82 MeV, fo, = 233 ± 49 MeV for (pole/pole), 

fot = 346 ± 82 MeV, /d, = 255 ± 54 MeV for (pole/const.). (50) 

From (^Tj), ( PI ) and ( |i5| ) the ratio of the vector and pseudoscalar decay constants fo*/ fn^ 
is given by 



1.36 



\Vc,Fr{ml)\ 



BD, 



m 



'6(3° 



D+D-*) 



B{BO D+Dj) 



1/2 



0.87 
096 



(51) 



which gives 



1.56 ±0.49 for (pole/pole), 
1.43 ± 0.45 for (pole/const.) 



(52) 



The decay constant is changed according to the behaviour of the form factor Fo{q^). 
However the amount of change is less than 10% as shown in (0). From this we know that the 
decay constant is not so much dependent on the behaviour of the form factor. Also when we 
consider the uncertainty from non-factorizable effects, the decay constant is changed within 10% 
discrepancy. In table H we show the results of fn* , and fo* / fos for different non-factorizable 
contributions. 

As discussed in , when we consider the penguin contributions with non-factorizable effects, 
the value of the decay constant is increased by 8%, however, for //^^ it is increased by up 
to 19%. So the ratio fo*/ foa is decreased by 9%. In table 4 we summarized the values of 
decay constant /z)., fr,^ and the ratio of foil fos from various sources. Our result for /^j^ 
agrees well with other theoretical calculations and experimental results within errors. For 



the ratio fD-ffos^ our results have a value greater than 1, however, Browder et al. has 
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a value less than 1. It seems that this ratio is more likely to be greater than 1 when we 
consider that the decay constant of p meson is 1.5 times greater than that of vr meson. The 
difference of the results by Cheng and Yang comes from the different method and using 
different Wilson coefficients. Their values come by comparing two non-leptonic decay modes, 
for instance B{B DDs{D^*W*))/B{B Dn). 



5. Conclusion 

By including the penguin contributions and the non-factorizable term contributions, we 
calculated the weak decay constants and fo* from D^i^u and B^ D^D^*\ In 

our analysis, we consider the QCD-penguin and Electroweak-penguin contributions in hadronic 
two body decays within the NDR renormalization scheme at next-to-leading order calculation. 
We also considered the effect of two different g^-dependence of the form factor for F^^(g^). 
The value of fu^ is changed by less than 10% for different form factors. 

The penguin effects for B DDg decay is quite sizable, and we obtained /d^ = 233 ± 49 
MeV for the monopole type of F^^, fo^ = 255 ± 54 MeV for the constant Ff^^ . When we 
considered the non-factorizable contributions, we obtained fo* = 346 ± 82 MeV for the D* 
meson. These values will be improved vastly when the large accumulated data samples are 
available at the Belle and BaBar experiments in near future. 
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B{B^ D+p-) 


B{B° D+K-*) 




B{B^ D+K~) 




xlO^ 


xlO^ 


xlO^ 


xlO^ 


(pole/pole) 


9.01 ± 1.54 


4.62 ± 0.79 


3.58 ±0.61 


2.74 ±0.47 


(pole/const.) 


9.01 ± 1.54 


4.62 ± 0.79 


3.57 ±0.61 


2.71 ± 0.46 


Experiments 


8.4 ± 1.6 ±0.7 




3.1 ±0.4 ±0.2 





Table 1: The obtained values of the branching ratios with oi = 1.02 and experimental mea- 
surements. 



Coefficients 


Real Part 


Imaginary Part 




1.168 


0.0 




-0.365 


0.0 


cf' 


2.25 -10-2 


4.5 -10-=^ 




-4.58 -10-2 


-1.36 -10-2 




1.33 -10-2 


4.5 ■lO--'^ 




-4.80 -10-2 


-1.36 -10-2 




2.37 -10-^ 


-2.88 -10-^ 




4.30 -10-^ 


0.0 




-1.11 -10-2 


-2.88 -10-^ 




3.75 -10-3 


0.0 



Table 2: The values of the effective wilson coefficient C^-'-' with the fi — mb{mb) = 4, 3 GeV, 
TTici^b) ~ 0.95 GeV in the NDR scheme at NLO calculation. 
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Toll = 2. mLU = 2 


Toll = 2. TOi.« = 5 


TOi.i. = 3. TOLi^ = 3 


Coeffs. 


Real Part 


Imag. Part 


Real Part 


Imag. Part 


Real Part 


Imag. Part 


ai 


0.985 


0.0 


0.985 


0.0 


1.046 


0.0 




0.219 


0.0 


0.219 


0.0 


0.024 


0.0 


03 


-4.00-10-^ 


-2.30 -10-3 


-4.00-10-^ 


-2.30 -10-3 


7.23-10-3 


-3.30 -10-5 


04 


-3.46 -10-2 


-1.14-10-2 


-3.46 -10-2 


-1.14-10-2 


-3.83 -10-2 


-1.12-10-2 


05 


-1.07 -10-2 


2.3 -10-=^ 


3.70 -10-=^ 


1.78 -10-3 


-2.70 -10-3 


-3.33 -10-5 


ae 


-4.13 -10-2 


-1.14 -10-2 


-4.53 -10-2 


-1.27 -10-2 


-4.36 -10-2 


-1.21 -10-2 


a-j 


-2.19 -10"^ 


-2.88 -10"^ 


-1.51 -10"^ 


-2.88 -10-^ 


-9.35 -10-5 


-2.88 -10-^ 


as 


3.11 -10-^ 


-1.44 -10-^ 


-3.82 -10-^ 


-5.77 -10-^ 


3.51 -10-^ 


-9.61 -10-5 


ag 


-9.27 -10-=^ 


-2.88 -10-^ 


-9.27 -10-=^ 


-2.88 -10-^ 


-9.90 -10-3 


-2.88 -10-^ 




-1.82 -10^^ 


-1.44 -10"^ 


-1.82 -lO^^ 


-1.44 -10"^ 


3.39 -10-5 


-9.61 -10-5 



Table 3: The values of the effective coefficients with /i — mb{mh) = 4, 3 GeV and mdm})) = 
0.95 GeV in the NDR scheme at NLO calculation. and a2i_i are defined by 021-1 = C2i-i + 
Cl(^ /Nt^^ and a2i = C^l^ + Cl{lj N^^f . Here we have taken {N^)ll for (V-A)(V-A) interaction 
and {Nc)lr for (V-A)(V-FA) interaction. 
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cILL 



c LR 



foi (MeV) 



fo. (MeV) 




346 ± 82 



231 ±48 



346 ± 82 



233 ± 49 



325 ± 77 



216 ±45 



1.57 ±0.50 



1.56 ±0.49 



1.57 ±0.50 




346 ± 82 



252 ± 53 



346 ± 82 



255 ± 54 



325 ± 77 



235 ± 50 



1.44 ±0.45 



1.43 ± 0.45 



1.44 ±0.45 



Browder et a/.[|T7| 
Hwang and Kim []TB| 
Cheng and Yang [|l^ 



243 ± 70 
362 ± 15 
266 ± 62 



277 ± 77 
309 ± 15 
261 ±46 



0.88 ±0.35 
1.17±0.02 
1.02 ±0.30 



Capstick and Godfrey pO| 
Dominguez |22 



UKQCD ^ 



BLs m 



MiLc m 



Becirevic et al. [E^ 



272 ± 16t!^, 



20 



290 ± 20 
222 ± 48 

230±7±35 
199 ± Stfltl' 
231 ± Utl 



1.18 ±0.18 



WA75 
CLEO 1 13 



CLEO 2 |2g 
BES |29| 
E653 [|0| 



238 ± 47 ± 21 ± 48 
282 ± 30 ± 43 ± 34 
280 ± 19 ± 28 ± 34 

mtlfo ± 40 

190 ± 34 ± 20 ± 26 



Table 4: The obtained values of fo* (MeV) and (MeV), and their ratio fo*/ fos^ the 
results from other theoretical calculations and existing experimental results. Here we refered 
the corrected fo^ values |^ for the experimental data - [^. 
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